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Abstract 

Context. Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation 

studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. 

Aims. We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical 

prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho H2O. 

Methods. In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 

MHz, corresponding to a line-of-sight velocity of 0.06 km s , which we employed to obtain the strength and linear polarisation of maser spikes 

in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut fiir 

Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. 

Results. We report the first astronomical detection to date of H2O maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 

GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected 

from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised 

maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the 

core of the line, they rise up to approximately 6% in its wings. 

Key words, stars: AGB and post-AGB - stars: winds, outflows - supergiants - circumstellar matter - masers - submillimetre 
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1. Introduction 

VY Canis Majoris is a highly luminous, strongly obscured, vari- 
able supergiant with a high infrared excess. Only ~ 1 % of the 
total luminosity is observed at optical wavelengths. The star's 
distan ce has been measured to be D — 1.1 kpc (Choi et al., 
2008) implying a luminosity L = 3 x 1O 5 L (Menten et al., 



2008). Three thermally emitted mid-infrared water vapour emis- 
sion lines seen in spectra obtained with the Short Wavelength 
Spectrometer (SWS) on the Infrared Space Observatory (ISO) 
indicate a mean radial velocity of order 20 ±2 km s _1 and a 
25 km s H2O outflow velocity (Neufeld et al., 1999), si gnific- 
antly l ower than the 32 km s _1 velocities that lReid & Muhlemanl 
J 19781) reported for 1612 MHz OH maser outflow. 

The effective tem perature of the central star, T* = 2,800K 
(Monnier et all fl99 9). combined with distance and inferred lu- 
minosi ty fixes the star's r adius at R* ~ 10 AU. Working at 
1 1 jUm. lDanchi et al.l (1 19941) reported a variable photospheric ra- 
dius ranging from 9.5 to 1 1 mas, where 10 mas corresponds to a 
stellar radius of 10 AU. 

VY CMa is a strong source of 616 — 523 22.235 GHz wa- 
ter vapour maser emission. Recent observations at the Atacama 
Pathfinder Experiment (APEX) telescope have revealed H2O 
maser emission also at eight submillimetre frequencies, ranging 
from 321 to 475 GHz (Menten et all 2008). Although a the- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



oretical model of Neufeld & Mel nickl d!991l) predicted that the 
22.235 GHz transition sh ould be far more lum inous than these 
submillimetre transitions, iMenten et alj d2008l) find them to be 
comparable in flux density, or at most a factor of <~ 6 lower. This 
suggests a need for further observations that might explain the 
discrepancies. 

One of strongest su bmillimetre masers predicted by 
iNeufeld & Melnickl I (1199 ll) is due to the 5 32 -4 4 i 620.701 GHz 
transition of ortho water, whose photon luminosity was expected 
to be roughly 16% that of the 22.235 GHz maser. Given the mul- 
tiple masers observed at 22.235 GHz we expected the 620.701 
GHz maser of VY CMa to exhibit similar multiple peaks at 
photon densities ~ 16% of those exhibited at 22.235 GHz. 



2. HIFI Observations at 620.701 GHz 

The Herschel Heterodyne Instrument for the Far Infrared (HIFI) 
cover s seven frequency band s, ranging from 488.1 to 1901.8 
GHz (IdeGraauwet all l2010t) . The 620.701 GHz radiation of 
the 532 — 441 maser of ortho H2O is observed in HIFI's Band 
IB, which nominally covers the range from 562.6 to 628.4 GHz. 
Two channels cover each of the frequency bands, one sensitive 
to linearly polarised radiation roughly parallel to the spacecraft 
horizontal (H) direction, the other roughly parallel to the vertical 
direction (V). For Band IB the H direction of polarisation is at an 
angle of 82.5° relative to the spacecraft V axis, while the V dir- 
ection of polarisation is at an angle of -7.5° to that axis. Our data 
were obtained in the HIFI high-resolution spectroscopy (HRS) 
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mode with spectral resolution 0.125 MHz, or line-of-sight velo- 
city resolution 0.06 km s _1 . 

The H and V beams on HIFI are not fully coincident. In 
Band IB they are separated by ~ 19% of the ~ 34.4" full-width- 
half-power beam diameter; the offset between the two beams is 
~ 6.6". V Y CMa is a spatia lly unresolved source at these fre- 
quencies (Dec in et al.L 12006). In our observations the star was 
positioned half-way between beam centres, i.e., displaced ~ 3.3" 
relative to the centre of each beam, although a random ~ 1.8" 
pointing error can slightly increase alignment uncertainty. 

In order to determine the orientation on the sky of any ob- 
served linear polarisation, a source has to be viewed at least 
at two different rotation angle s relative to the tele scope. The 
Herschel Space Observatory (Pilbra tt et al.L |2010), however, 
cannot be rotated without producing undesirable thermal drifts. 
A rotation is best achieved by observing a target at two epochs 
separated by a number of weeks. The further the target lies above 
the ecliptic plane, the faster is the rotation produced. In our ob- 
servations, an interval of three weeks between two sightings of 
VY CMa resulted in a rotation of the telescope of ~ 16° about 
the line of sight to the star. 

We first observed the star on March 21, 2010 in two con- 
tiguous segments, lasting 10,139 seconds apiece, for a total 
observing time of 5.633 hours. The first segment began at 
08:48:38.0 UT; the second was started at 1 1:45:00 UT, 7 minutes 
and 23 seconds after the first had terminated. Our second epoch 
of observations started at 17:02:33 UT on April 11, 2010 and 
again ran for 10,139 seconds or 2.816 hours. The system tem- 
perature for these observations ranged between approximately 
80 K and 100 K with respective RMS noise levels for the first 
and second sets of observations of 6. 1 and 5.3 Jy for the H chan- 
nel, and 4.3 and 4.5 Jy for the V channel at the resolution of 
0.125 MHz. 



3. Polarisation Analysis of the 620.701 GHz Maser 
Components 

Linear polarisation measurements are often optimised to determ- 
ine th e Stokes Q and U parameters indepen dently from each 
other dLi et all 120081; IHezareh & Houdd.120101) . If we define four 
intensity measurements I<p within the equatorial system with q> 
set to 0° when pointing north and increasing eastwards, then 



These 2/V intensities yield the polarisation if we first define 



/ = (/o+/ 45 +/90 + /l35)/2 

Q = I0-I90 
u = /45-/135. 



(i) 

(2) 
(3) 



Although the state of linear polarisation is completely defined 
with these equations, it is also commonly expressed with 
the polarisation fraction p = \/Q 2 + U 2 /I and angle 6 = 
0.5arctan(£//(2). Alternatively, we can also write 



Q = pi cos (29) 
U = pi sin (20). 



(4) 
(5) 



With HIFI we do not generally have access to four independ- 
ent measurements precisely sampled with a spacing of 45°, as in 
Equations dXJ-([3j. Instead, we may have access to N irregularly 
spaced observing epochs with the HIFI vertical polarisation axis 
(V) oriented at an angle <p, relative to north (i = 1,2,..., TV) and 
the horizontal polarisation axis (H) at (f>, + 90°, each observation 
yielding a pair of intensities Iy i and /#. . 



dt = Iy i - 1 Hi 

= Qcos(2(pi) + Usm{2(pi). 



(6) 



Introducing a vector d with elements di/ot, where d\ is the in- 
tensity difference of Equation d6) and <7,- its uncertainty (i.e., the 
corresponding noise intensity), we can write the matrix equation 



d = A s, 

where the matrix A and vector s are given by 
cos(2<pi)/(Ti sin(2<pi)/o"i 

cos (2<pjv) /ojv sin (2<p#) /on 

Q 

u 



(7) 



(8) 



(9) 



The vector s, and therefore the Stokes parameters, are eas- 
ily obtained by inverting Equation Q using the pseudo-inverse 
matrix of A. We then find 



(A r A)~V 



d. 



(10) 



The elements of the covariance matrix C = (A r A) can be 
shown to yield the corresponding uncertainties on the estimated 
Stokes parameters with 






(ID 



where o2 = (sfSj) — (${) (sj) with i,j = QorU and Si the ap- 
propriate component of the vector s. It is also possible to determ- 
ine the uncertainties o p and Og for the polarisation fraction and 
angle, respectively, with 



-^ (Q 2 a 2 Q + U 2 a 2 v + 2QUa 2 QU 



4p 4 I 



-AQ 2 <yhv 



-U 2 o 2 QQ -2QUc; 2 QU ) 



(12) 
(13) 



We have applied this polarisation analysis to our observa- 
tions of the 620.701 GHz H 2 5 32 -4 41 ortho-transition. Our 
results shown in Figure[T]were obtained from the two aforemen- 
tioned observations made at epochs of corresponding position 
angles of 261.27° and 277.46°. 

Although there are no obvious strong polarisation signals 
from the maser emission peaks, we clearly detect polarisa- 
tion levels ranging from p ~ 1.5% to p ~ 6% in regions of 
significant line intensity (i.e., from approximately -5 to 45 
km s _1 ). Furthermore, the observed anti-correlation of the po- 
larisation fraction with the Stokes I intensity is sim ilar to 
previous ground-based po larisation observations dGirart et all 
120041: IHezareh & HoudeL l2010h aimed at the detection of 
the Goldreich-Kylafis effect in non-masing molecular lines 
(iGoldreich & Kvlafisl Il98lt ICortes et all 120051). which appears 
to have first been detected in evolved stars (iGlennet a l.. 1997). 
We will discuss the relevance of the Goldreich-Kylafis effect for 
our observations in section 7 below. 
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Figure 1. Polarisation spectrum of the 620.701 GHz H2O 532 — 
441 ortho-transition for VY-CMa. Shown are, from top to bot- 
tom, the Stokes / and polarisation fraction (symbols and right 
axis), the polarization angle, and the (unsigned) polarised flux. 
The polarisation fraction and flux were corrected for posit- 
ive bias but not for instrumental polarisation. A simple linear 
baseline was removed from each V and H channel observation 
and all polarisation vectors shown satisfy p > 3a p and / > 10(7/ 
at the corresponding velocities, where (7/ is the uncertainty (or 
noise level) on the Stokes / spectrum (i.e., the average of the 
h +Ih spectra measured at the two epochs). The polarisation 
fraction and angle have their spectral resolution reduced by a 
factor of 20 compared to Stokes /, which is at the maximum res- 
olution. 



4. Results 

Several instrumental capabilities of HIFI and the observations 
obtained with them may be noted: 

(i) The two orthogonally polarised HIFI receivers are well 
matched and extremely stable. But observations of extended 
sources need to be conducted with caution. A slight misalign- 
ment of H and V receivers can lead to a "false polarisation" that 
reverses polarity at half-year intervals (see Appendix lAl). 

(ii) The misalignment of the HIFI receivers does not appear 
to affect observations of unresolved sources. Our observations 
realised at two observing epochs indicate that instrumental po- 
larisation, which could be in part due to errors in the relative 
calibration between the two receiver chains of Band IB, cannot 
exceed a measure of order 1 — 2% (see section 7 below). 




Figure 2. Residual spectrum resulting from the subtraction of 
the 620.701 GHz Stokes / spectra taken 3 weeks apart. The re- 
sidual signal has an amplitude < 1%, underlining the accuracy 
of the relative calibrations as well as the remarkable stability of 
the masers. 



(iii) The polarisation of VY CMa is not significant near 
the peak of the 620.701 GHz H2O 532 —4-a\ line, but rises up 
to ~ 6% in the wings of the spectrum in a manner consistent 
with polarisation due to the Goldreich-Kylafis effect discussed 
in greater detail in section 7 below. 

(iv) The stability of the 620.701 GHz masers is remarkable. 
The variation over a three week period is <^ 1% (see Figure 2). 

(v) As Figure [3] shows, the spectral profile of the 620.701 
GHz and 22.235 GHz masers appears remarkably similar, ex- 
cept that the relative expansion velocities between peaks, along 
the line of sight, is a factor of 2.3 greater at 620.701 than for con- 
temporaneous observations at 22.235 GHz (see AppendixlBl). 

(vi) Assuming the 620.701 GHz masers to sit atop a broad 
pedestal with flux density ~ 400 Jy, we find the main 620.701 
GHz maser peak flux density to be ~ 800 Jy, compared to ~ 
1900 Jy at 22.235 GHz. The velocity spread of the expanding 
620.701 GHz maser peaks is ~ 2.3 times wider than at 22.235 
GHz, implying a 620.701 GHz photon luminosity roughly 2.4 
times lower than at 22.235 GHz — comparable in luminosity to 
the submillimetre masers observed by Menten et al. (2008), and 
about 2.5 times more luminous than lNeufeld & Mel nickl (Il99lh 
predicted. 



5. Discussion 

A spherical ly symmetric mod el for the outflows of VY CMa de- 
veloped by iDecin et alj (2006) envisions dust formation in the 
cooling outflow from the star at distances ~ 10/?*. Radiation 
pressure then accelerates the dust, and with it also the ambi- 



ent gas, from velocities <J 5 km s to ~ 25 km s at a dis- 
tance of 20/?*. The kinetic temperature at 20/?* is of the or- 
der of - 1000 K, sufficient to excite both 620.701 and 22.235 
GHz masers. We propose that the masing outflows we observe 
propagate along several mutually common directions from the 
star, and that the faster expanding 620.701 GHz masers lie at 
greater distances from the star than the central cluster of masers 
peaking at 22.235 GHz. We picture this outflow consisting of gas 
ejected obliquely to the line of sight to VY CMa but fanning out 
over a considerable angular width to account for the observed 
range of outflow velocities. Common ejection of both the denser 
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Figure 3. Comparison of the VY CMa 620.701 and 22.235 GHz 
masers, both observed on April 11, 2010. The 620.701 GHz 
(top) and the 22.235 GHz (bottom) spectra have different velo- 
city scales, and their peaks fall at slightly different vi sr velocities 
that have been shifted to coincide. Superposed on the 620.701 
GHz spectrum is a Gaussian pedestal subtracted in the spectrum 
just below to display the similarities of the 620.701 and 22.235 
GHz maser spectra. Not knowing the source or true shape of the 
pedestal, we kept its spectrum simple. 



masing regions and the less-dense ambient medium responsible 
for the pedestal emission can then account for the observed velo- 
city ranges of both these outflow components. The region from 
which the 22.235 GHz radiation originates may be at the inner 
radius of the flow; for, Neufeld & Melnick (1991) find that, for 
comparable gas densities and H2O abundances, maser emission 
at 22.235 GHz should be maximized at line-of-sight velocity 
gradients ~ 30 times lower than those maximizing 620.701 GHz 
emission. 

A comparison of our 620.701 GHz spect rum to the spectra of 
the star's submillimeter masers observed by Me nten et al.l (12008) 
shows that their spectra, observed in June 2006, four years earlier 
than ours, generally exhibited only two peaks separated by <; 6 
km s _1 , whereas our spectrum at 620.701 GHz clearly exhib- 
its four peaks ranging over separations of order ~ 25 km s . 
However, their 22.235 GHz spectrum also looks quite differen t 
from ours. This is not surprising since Esimbek & Zheng! (12001) 
found that, between August 1993 and August 1999, the shape of 
the 22.235 GHz spectrum significantly changed and peak emis- 
sion shifted from vi sr ~ 17 to 34.2 km s . 



Nevert heless, a feature common to all the submillimetre 
spectra of lMentenet alJ (120081) is a broad pedestal that dimin- 
ishes in strength with rising upper-excitation temperature E u of 
the observed transitions. For E u ~ 725 K it remains strong with 
wings extending over a full width > 40 km s , whereas at E u 
= 1861 K, the pedestal has entirely disappeared. The radiation 
emitted in the pedestal thus appears to be due to collisional ex- 
citation by gas flowing out from the star at spee ds > 20 km 
s -1 cl ose to the terminal velocity in the model of iDecin et al.l 
(120061) . The pedestal we observe for the transition at 620.701 
GHz, whose E„ = 732 K, has the same broad wings extending 
roughly symmetrically to either side of the maser outflow peaks. 

Low E u levels thus excited can then radiate in three ways, 
spontaneo usly, through stimulated emission, or induced by colli- 
sions. As Neufeld & Melnick ( 1991) point out, stimulated emis- 
sion leads to masing only when the ratio of water vapor density 
to the velocity gradient along the line of sight is sufficiently high. 
We interpret the maser peaks in our 620.701 GHz spectrum as 
emanating from denser clumps in the VY CMa outflow within 
which the velocity gradient is low, and the 620.701 GHz pedes- 
tal as emission from an ambient medium characterized by lower 
density and/or higher velocity gradients, and thus unable to ra- 
diate significantly through stimulated emission to sustain maser 
amplification. This ambient medium, and its corresponding radi- 
ation, will bear the imprint of the Goldreich-Kylafis (GK) effect 
when also subjected to anisotropic radiation or optical depths, 
and permeated by a weak magnetic field. 

Although the low polarization levels we detect are expec- 
ted for water masers (Watson, 2009), the rise in polarisation 
seen in the wings of the spectrum in the top panel of Figure 
Q] deserves careful consideration. If we assume that the masers 
simply maser-amplify the seed radiation emanating from the am- 
bient gas responsible for the pedestal component, then we may 
have a qualitative explanation for the shape of the polarisation 
spectrum. That is, the polarization fraction should be approxim- 
ately the same at the position of, and between, maser peaks, as 
observed in Figure [TJ and should rise in the wings of the line 
where the optical depth of the pedestal component is lowest, 
in a manner consistent with the GK effect (Cortes et al., 2005). 
Although maser-amplification of GK polarisation could lead to 
different levels of polarisation at the maser peaks, other factors 
(e.g., magnetic field orientation) could also suppress this effect. 
We also note that in cases where the pedestal component was 
highly saturated, the level of polarisation measured in its core 
woul d provide a measure of th e amount of instrumental polarisa- 
tion dHezareh & Houdel 12010). Although this scenario probably 
does not perfectly fit our observations, the polarisation fraction 
measured in the vicinity of the line core suggests that the instru- 
mental polarisation is approximately 1 - 2%. 



6. Conclusions 

We report the first astronomical observations of the 620.701 GHz 
532—441 submillimetre maser of ortho^O. The maser peaks in 
VY CMa show weak linear polarisation at levels consistent with 
the Goldreich-Kylafis effect across the spectral line profile. 

Acknowledgements. HIFI has been designed and built by a consortium of 
institutes and university departments from across Europe, Canada and the 
United States under the leadership of SRON Netherlands Institute for Space 
Research, Groningen, The Netherlands and with major contributions from 
Germany, France and the US. Consortium members are: Canada: CSA, 
U.Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, 
MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio 
Astrofisico di Arcetri- INAF; Netherlands: SRON, TUD; Poland: CAMK, 



Martin Harwit et al.: Polarisation Observations of VY Canis Majoris H2O 532 — 441 620.701 GHz Maser Emission with HIFI 



CBK; Spain: Observatorio Astronmico Nacional (IGN), Centro de Astrobiologia 
(CSIC-INTA). Sweden: Chalmers University of Technology - MC2, RSS & 
GARD; Onsala Space Observatory; Swedish National Space Board, Stockholm 
University - Stockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: 
Caltech, JPL, NHSC. Support for this work was provided by NASA through an 
award issued by JPL/Caltech. We thank the HIFISTARS consortium for permis- 
sion to use their VY CMa 556.9 GHz data for calibration purposes. We would 
like to acknowledge that Tom Phillips first pointed out to us that HIFI might 
yield useful polarisation observations. And, finally, we extend thanks to the an- 
onymous referee for perceptive insights and suggestions that greatly improved 
this paper. 



References 

Attai'd, M., Houde, M., Novak, G., and Vaillancourt, J. E., 2008, PASP, 120, 805 
Bachiller, R., Perez Gutierrez, M., Kumar, M. S. N., and Tafalla, M., 2001, A&A, 

372, 899 
Choi, Y. K., et al., 2008 . PASJ, 20, 1007 

Cortes, P. C, Crutcher, R. M., and Watson, W. D. 2005, ApJ, 628, 780 
Danchi, W. C, Bester, M., Degiacomi, C. G, Greenhill, L. J., and Townes, C. H. 

1994, AJ, 107, 1469 
Decin, L., Hony, S., Justtanont, K., Tielens, A. G G. M., & Waters, L. B. F. M., 

2006, A&A, 456, 549 
de Graauw, Th. et al., A&A, in press 

Esimbek, J. & Zheng, X.-W., 2001, Chinese Physics Lett., 18, 458 
Girart, J. M., Greaves, J. S., Crutcher, R. M., and Lai, S.-P. 2004, Ap&SS, 292, 

119 
Goldreich, P., and Kylafis, N. 1981, ApJ, 243, L75 

Glenn, J., Walker, C. K., Bieging, J. H, and Jewell, P. R. 1997, ApJ, 487, L89 
Hezareh, T., and Houde, M. 2010, PASP, in press (arxiv.org/abs/1005.0881) 
Li, H., Dowell, C. D., Kirby, L., Novak, G, and Vaillancourt, J. E. 2008, 

Appl. Opt., 47, 422 
Menten, K. M., Lundgren, A., Belloche, A., Thorwirth, S., and Reid, M. J., 2008, 

A&A, 477, 185 
Monnier, J. D., Geballe, T. R., and Danchi, W. C. 1999, ApJ, 512, 351 
Neufeld, D. A., and Melnick, G. J. 1991, ApJ, 368, 215 
Neufeld, D. A., Feuchtgraber, H, Harwit, M., and Melnick, G. 1999, ApJ, 517, 

L147 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., Crone, G., Doyle, D., Gageur, U, 

Heras, A. M., Jewell, C, Metcalfe, L., Ott, S., and Schmidt, M. 2010, A&A, 

in press 
Reid, M. J., and Muhleman, D. O. 1978, ApJ, 220, 229 
Vlemmings, W. H. T., Diamond, P. J. & van Langevelde, H. J. 2002 A&A, 394, 

589 
Watson, W. D. 2009, RevMexAA (serie de Conferencias), 36, 113 

1 Cornell University, Center for Radiophysics & Space Research, 
51 1 H street, SW, Washington, DC 20024-2725 USA 
e-mail: harwit@verizon.net 

University of Western Ontario, Department of Physics and 
Astronomy, London, Ontario, Canada N6A 3K7 

Johns Hopkins University, Department of Physics and Astronomy, 
Baltimore, MD, 21218 USA 

4 IPAC, Caltech, Pasadena, CA 91925 USA 

5 Consejo Superior de Investigaciones Cientificas, 28006, Madrid, 
Spain 

6 Katolieke Universiteit Leuven, Institut voor Sterrenkunde, 
Heverlee 3001 Belgium 

7 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, 
53121 Bonn, Germany 

8 National University of Ireland, Maynooth, Department of 
Experimental Physics, County Kildare, Ireland 

9 Space Research Organization of the Netherlands (SRON) 9700 
AV Groningen, Netherlands 

10 University of Waterloo, Department of Physics and Astronomy, 
Waterloo, Ontario, Canada N2L 3G1 

11 Harvard-Smithsonian Center for Astrophysics Cambridge, MA 
02138 USA 

12 European Space Astronomy Centre, Urb. Villafranca del Castillo, 
P.O. Box 50727, Madrid 28080 Spain 

13 SOFIA Science Center, Universities Space Research Association, 
NASA Ames Research Center, Moffett Field, CA 94035-0001 USA 

14 Observatorio Astronomico Nacional (IGN), Alfonso XII N°3, E- 
28014 Madrid, Spain 



15 Observatorio Astronomico Nacional (IGN), Ap 112, E-28803 
Alcala de Henares, Spain 

16 Sterrenkundig Instituut Anton Pannekoek, University of 
Amsterdam, Science Park 904, NL-1098 Amsterdam, The 
Netherlands 

17 Department of Astrophysics/IMAPP, Radboud University 
Nijmegen, Nijmegen, The Netherlands 

18 Onsala Space Observatory, Dept. of Radio and Space Science, 
Chalmers University of Technology, SE-43992 Onsala, Sweden 

19 The Netherlands and Astronomical Institute, Utrecht University, 
Princetonplein 5, 3584 CC Utrecht, The Netherlands 

20 Department of Astronomy, AlbaNova University Center, 
Stockholm University, SE-10691 Stockholm, Sweden 

21 Joint ALMA Observatory, El Golf 40, Las Condes, Santiago, 
Chile 

N. Copernicus Astronomical Center, Rabianska 8, 87-100 Toruri, 
Poland 



Appendix A: A Cautionary Note on Observations 
of Extended Sources 

As part of the Performance Verification of HIFI in space, spec- 
troscopic observations of the 12 C 16 7 = 5^4 576.3 GHz 
transition were obtained at the position of the chemically-active 
bow shock Bl drive n by the LDN 1157 Class protostar 
(lBachilleretal.Ll200lh . Initial observations in HIFI's Band IB 
were made on August 1, 2009; a second set was obtained 186 
days later on February 4, 2010. The spectral shape of the source 
remained essentially unchanged but initially the signal from the 
high-frequency edge of the line appeared stronger in the V direc- 
tion, when the H and V peaks were matched, whereas half a year 
later it appeared stronger by the same ratio in the H direction. 

This "false polarisation" effect was due to the small mis- 
alignment of the H and V beams described earlier, which im- 
aged slightly different portions of a scene onto the two receivers 
( Attard et al., 2008). Unless otherwise specified, observations of 
a targeted source are centred on a position halfway between the 
two beam centres. Because the viewing direction of Herschel 
at all times is constrained by the need to keep the plane of the 
sun-shield roughly perpendicular to the radius vector to the Sun, 
sources close to the ecliptic plane can only be viewed at half year 
intervals. Over this interval the telescope aperture rotates 180° 
about its viewing direction on the sky. Thus, portions of LDN 
1 157 Bl initially viewed by the V channel were viewed, half a 
year later, by the H channel and vice versa. It is clear that this 
effect was not due to linear polarisation because a 180° rotation 
of the telescope leaves the polarisation direction unchanged. 

For VY CMa, where the entire source was relatively well 
centred on each beam, this effect was negligible. But for linear 
polarisation observations of an extended source the best strategy 
may be to make a small map with HIFI over the area of interest, 
where the spectra obtained in the H and V channels at each pos- 
ition of the source, can be individually and directly compared. 

Appendix B: Ground-based Observations at 22.235 
GHz 

Overlapping with the second epoch of 620.701 GHz observa- 
tions, we observed the VY CMa 22.235 GHz maser with the 
Effelsberg 100-meter telescope of the Max-Planck-Institut fiir 
Radioastronomie. Those observations were begun on April 11, 
2010 at 16:40 UT and lasted until 17:20 UT and are shown in 
Figure[3] Data in two linear polarisations were obtained with the 
two channels of the K-band or 1.3 cm receiver at the primary 
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focus. The frequency resolution was 6.104 kHz, corresponding 
to a velocity resolution of 0.082 km s _1 . Although these ob- 
servations were realised at only one epoch and are, therefore, 
insufficient for a full characterisation of t he linear polarisation 
state, a search by Vlemmings et alj d2002l) found that the linear 
polarisation of the VY CMa 22.235 GHz masers is well below 
1%. 



